ITP-SB-95-19 
May, 1995 



The soft pomeron and nonperturbative corrections 
in quark-quark scattering. 



Irina A. Korchemskaya * 

Institute for Theoretical Physics, 
State University of New York at Stony Brook, 
Stony Brook, New York 11794 - 3840, U.S.A. 

and 

Universitd di Parma and INFN, 
Gruppo Collegato di Parma, 1-43100 Parma, Italy 



Abstract 

We consider elastic quark-quark scattering at high energy with fixed momentum trans- 
fer t and perform factorization of soft-gluon exchanges into a vacuum expectation value 
of Wilson lines. Taking into account nonperturbative corrections whose structure is pre- 
dicted from infrared renormalon analysis, we represent the scattering amplitude as an 
asymptotic series. In the region of small momentum transfer, where the nonperturbative 
corrections are dominant, the scattering amplitude is Gaussian distribution in t with a 
slope depending on a nonperturbative scale. A nonperturbative origin of the soft pomeron 
is thus identified. 
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1. Introduction 



Regge theory [1] explains a large class of experimental results for hadronic scattering amplitudes 
at high energy s and fixed transferred momentum t. After almost 35 years, however, it remains 
a challenge to understand Regge theory within the framework of fundamental quantum field 
theory. One such attempt led to the development of the BFKL pomeron [2], which describes 
the compound state of two rcggcizcd gluons with vacuum quantum numbers. Found in the 
leading logarithmic approximation, the BKFL pomeron, sometimes called the hard pomeron, 
leads to scattering amplitudes which violate the Froissart bound. Recently [4], [5], [6] it was 
shown that in generalized leading logarithmic approximation QCD is described by an effective 
2-dimensional field theory, which is equivalent to an XXX Heisenberg magnet of spin s — 0. 
This theory takes into account the propagation of infinite numbers of interacting Reggeons in 
the t-channel and restores the unitary of the S-matrix. 

This paper is devoted to the "soft" pomeron [7] , [8]. While this idea is phenomenologically 
very successful [9], and we know from experiment the pomeron trajectory, a{t) = 1.08 + 
0.25t, a deep understanding how the pomeron appears in QCD is still lacking. It is widely 
believed that the soft pomeron has a nonperturbative origin. It was proposed in [10] to include 
nonperturbative effects through a modification of the gluon propagator. Another ideas, in 
particular the application of the method of the stochastic vacuum models were discussed in [11]. 
In [12] the high energy interactions were described in terms of a two-dimensional sigma-model 
action. The approach developed here follows from a different perspective. We shall exploit the 
viewpoint that perturbation theory itself can predict the form of nonperturbative corrections, 
and show that ambiguities of the perturbative series caused by infrared renormalons allow one 
to identify the structure of nonperturbative corrections. In case of e^e~ annihilation, which 
admits the operator product expansion, the nonperturbative corrections can be parameterized 
by local vacuum condensates. Further analysis [14] revealed that the perturbative series is not 
Borel summable. The singularities of the Borel transform, which are called IR renormalons [13] 
imply that the physical quantity will be well defined if the ambiguity caused by IR renormalons 
is compensated by an ambiguity in the definition of local vacuum condensates. In [15] the idea 
of IR renormalons was generalized to hadronic processes (jet cross sections, inclusive Drell-Yan 
lepton-pair production) to which the operator product expansion is not applicable. It was 
found that nonperturbative corrections are parameterized by new parameters, associated with 
vacuum expectation values of nonlocal operators involving Wilson lines and the gluon field 
strength. The strength of the nonperturbative corrections is determined by the position of the 
leading renormalon. The relation between infrared renormalons and power corrections was also 
discussed in [16], [17], [18], [19], [20]. 

Analyzing high-energy hadronic scattering, we consider hadrons as consisting of partons. 
The soft pomeron is thought to couple to the valence quarks only. This additive quark rule is 
supported by much experimental data. Quark-quark scattering naturally involves into consid- 
eration the Wilson lines. In nonperturbative QCD the Wilson line appears, for example, in the 
representation of the quark propagator S{x,y; A) as a sum over random paths between points 
x and y [22], [23] and it takes into account the interactions of gauge field with the color current 
created by the quark moving along the path. The remarkable fact, however, is that in high 
energy scattering, s ^ — t, the quarks move along the straight lines. In perturbative QCD the 
Wilson lines describes the infrared asymptotics of the quark propagator [21]. Moreover one can 
expand the Wilson line in powers of gauge fields and reproduce the eikonal approximation for 
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the interaction vertices of quark with soft gluons. 

Our strategy is the following. We start with factorization of the soft gluon exchanges into a 
vacuum expectation value of Wilson lines and represent the quark-quark scattering amplitude 
as an expectation value of a Fourier transformed Wilson line, evaluated along an integration 
path which consists of two semiclassical quark trajectories, separated by an impact parameter 
in the transverse direction [24]. For the case in which the quarks are described by hght-like 
Wilson lines the similar formula was proposed by Nachtmann [25] . We shall go on and find the 
expression for the quark-quark elastic scattering amplitude, which resembles an eikonal formula 
of Cheng and Wu [3]. We shall derive it, however, from the renormalization properties of the 
cross singularities of Wilson loops [24] . 

Calculated to the lowest order of perturbation theory the quark-quark scattering amplitude 
gets large perturbative corrections such as (ctslogslogi)". Resummation of these Sudakov 
corrections can be performed using an "evolution equation" technique [26], [27]. In the present 
paper, we generalize the method [24] to perform the resummation of both kind of corrections: 
perturbative and nonperturbative in the impact parameter space. The perturbative Sudakov 
corrections come from very soft virtual gluons with transverse momenta much smaller than 
the hadronic scale. In this region the nonperturbative corrections are very large. As we will 
show below, they have the form of power corrections (pAg^^^fe^)", where p is a parameter 
characterizing the nonperturbative interactions, while the perturbative corrections behave as 
(ttslogft^)". The dependence of the scattering amplitude on t comes from the dependence of 
the Wilson lines on the impact parameter b. This means that we expect a more or less realistic 
prediction for the t dependence of the scattering amplitude, and our soft pomeron predicts a 
linear Regge trajectory. It is problematic to find the intercept of the soft pomeron from our 
model and we can only say that it is close to unity. 

The paper is organized as follows. In Section 2 we review the application of the Wilson loop 
formalism to quark-quark scattering [24], [28]. In Section 3 we show that the Wilson loop in 
perturbation theory contains an ambiguous contribution to power corrections caused by infrared 
renormalons. In Section 4 using the prediction for the asymptotics of the cross anomalous 
dimension [28] in all orders of perturbation theory we include the nonperturbative corrections in 
the expression for the scattering amplitude, and represent the scattering amplitude as a Mellin 
integral. In Section 5 we calculate the scattering amplitude for the case of a frozen coupling 
constant. In Section 6 we consider the case of a running coupling constant and represent 
the scattering amplitude as an asymptotic series. Here we discuss the origin of the asymptotic 
expansion as a consequence of the presence of a singularity in the running coupling constant. We 
also derive the main results of this analysis: a Gaussian distribution over transferred momenta 
in the scattering amplitude, shrinkage of the distribution with increasing energy, and a crossover 
region in the differential cross-section. Section 7 contains concluding remarks. 



2. Scattering amplitude in perturbative QCD 

We consider near forward elastic quark-quark scattering at high energy and fixed transferred 
momentum in the following kinematics: 

s, > > > Aqcd ■ 

Here s = (pi+P2)^ is the invariant energy of quarks with mass m ,t = {pi—PiY is the transferred 
momentum and is an IR cutoff. Let us explain the origin of this kinematics. The quark- 
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quark scattering should be thought of as embedded in a physical process involving hadron- 
hadron elastic or inelastic scattering. Considered in isolation, the q-q scattering amplitude 
has IR divergences and an IR cutoff is necessary. In exclusive and inclusive physical processes 
involving the q-q scattering, the IR divergences are canceled and the IR cutoff is replaced by 
a dynamically generated transverse momentum scale of the hadronic state. Therefore for a 
purely perturbative calculation we should choose S> ^^qcd- To include nonperturbative 
effects, however, we find it useful to relax this condition such as > Aq^^^. As a particular 
example of hadronic state one can consider the perturbative onium state [29] built from heavy 
quarks with mass m. In this case the mass m has a meaning of the transverse size of the hadron. 

In the center of mass frame of the incoming quarks, quark momenta have the following 
light-cone components: pi = m/\/2(e'''/^, e"'''/^, 0) and p2 = m/\/2{e~'^^'^,e^^'^,Q). In the limit 
s ^ the angle 7 between quark velocities become large and both quark move close to 
the "+" and "-" light-cone directions. The components of the total momentum transfer k — 
Pi — p'l are = —k~ = 0{t/^/s) and k^ = 0{t). Thus, in the limit s ^ — t we can neglect 
the longitudinal components of transferred momentum and put k = (0'^,0~,A;). Since the 
transferred momentum is much smaller than the energies of incoming quarks, the quarks interact 
each other by exchanging soft gluons in the i-channel with total momentum k. Interacting with 
each of soft gluons quark does not alter its velocity in the limit —t <S and thus the only 
effect of its interaction is the appearance of an additional phase in the quark wave function. 
This phase, the so-called eikonal phase, is equal to a Wilson line P exp{i (ix^A^(a;)) evaluated 
along the classical trajectory C of quark in the direction of the quark velocity. We combine 
the eikonal phases of both quarks and obtain the representation for the scattering amplitude 
as [24]: 

V/i^, 5) = smh^Jcfbe-^-^Wff'i^,^)^) , t^-P, (2.1) 

where the line function Wj_j is given by 

Wi'/{-f,b'^X'^) = {0\T\Pexp(ig H davi ■ Aivia))]' [Pexp (i^ H dp V2 ■ A{v2p + b))]^ |0). 

L V J —00 /Jj L V J —00 ^ -I i 

(2.2) 

Here the line function Wl^-' contains color indices of both incoming (i, j) and outgoing («', j') 
quarks. The two Wilson hues are defined in the fundamental representation of the SU(N) gauge 
group and evaluated along infinite paths in the direction of the quark velocities vi — pi/m and 

— * 

1^2 = Vijfn. The integration paths are separated by impact vector b — (0+,0~,6) in the 
transverse direction, - b = V2 - b = ^. 

The scattering amplitude (2.1) depends on the quark velocities V\ and ?;2, the transferred 
momentum k and the IR cutoff A. These variables give rise to only two scalar dimensionless 
invariants: {y\V2) = s/m? and t/X^, as explicitly indicated in (2.1). The s-dependence of the 
amplitude comes from the dependence on the angle 7 between quark four- velocities Vi and V2, 
defined in Minkowski space-time as {V1V2) = cosh 7, while its t-dependence is related to the 
dependence of the line function on the impact vector b. In the limit of high-energy quark-quark 
scattering (s 3> m^) we have 

7 = log ^ > 1 . (2.3) 

The line function W^y is divergent for 6 = 0. This divergence, the so called cross divergence, 
has an ultraviolet origin, because for 6 = the integration paths of Wilson lines (See Fig. la) 
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cross each other. According to the general analysis in [30], the Wilson line (Wi)^/ of Fig. la 
is mixed under renormalization with the Wilson line {W2)ij of Fig. lb. As a consequence, the 
renormalized line functions Wi and W2 satisfy the following renormalization group equation: 

(/^l^ + Wa = -r«L(7, g)W, a, & = 1, 2, (2.4) 

where Wa = {Wa)ij and is the renormalization scale. Here, F^^^^g is the cross anomalous 
dimension, which is a gauge-invariant 2x2 matrix, depending only on the coupling constant 
and the angle 7 between the lines at the cross point. Solving the RG equation for Wi{'y, p-) 
with the boundary conditions W^i(7, 1) = Si'i6j>j and ^^2(7, 1) = SfiSiij and identifying the UV 
cutoff //^ with 1/6^ we find the following expression for the scattering amplitude [24], [28]: 

[-^^ 5) = / d%e-''-\An(% h^\^)5,i5^'j + Ai2(7, h^\^)5rA'i) , (2.5) 

where An and A12 are elements of the 2x2 matrix 

^(7, b'X') = Texp f^'^ Feross(7, as(r))^) . (2.6) 

We conclude that the asymptotic behavior of the scattering amplitude is governed by the matrix 
of the cross anomalous dimension rcross(7, Ois)- The expression (2.6) takes into account not only 
all logt and logs corrections but nonperturbative corrections as well. As follows from (2.6), 
the t-dependence originates from the evolution of the coupling constant. 
The one-loop expression for the matrix cross anomalous dimension is 

rc.oss(7,5) = -rcross(7), 
TT 

/ -fcoth7 mcoth7 \ 

icrossu; V-7coth7 + l + i7rcoth7 Ar(7coth7- 1) - f coth7y' • ^ > 

In the large s limit we have the following expressions for the eigenvalues of the matrix rcross(7)- 



e 2?7r S' 

N\og—-N- — + 0(log-i — ), (2.8) 



r- = ^'^^ log-^ 4 + o{iog-^ 4) . (2.9) 

The scattering amplitude can be decomposed into singlet and octet invariant amplitudes cor- 
responding to exchanges in the i-channel with quantum numbers of the vacuum and the gluon, 
respectively: 

Tff = 5,,5riT^ + i!t>it^,;T^. (2.10) 

Using the one-loop expression (2.7) for the matrix rcross(7,5') we find the following expressions 
for the invariant amplitudes: 

rW = -sinh7(^T+-T_), (2.11) 
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r 

where 



T^^) = 2i7rcosh7-^(T+ - T_) (2.12) 



In the leading logt and logs approximation the result for the invariant scattering amplitude 
has the standard reggeized form [31]: 

where Tsom = The nonleading logs corrections drastically change the behavior of 

the scattering amplitude. Indeed, the functions T± (2.13) have a Regge-like behavior, and 
for r+ ^ r_ we have T+ <^ T_. Thus, the high-energy behavior of the invariant scattering 
amplitudes (2.11), (2.12) is dominated by the contribution of T_, and as a consequence the 
amplitude of the octet exchange is suppressed by the factor l/r_|_ = O(log~^(s/m^)) compared 
to the amplitude of the singlet exchange. Therefore we will start from the main result of 
perturbative QCD: 

T(°)=sinh7r_ (2.15) 
T(8) = -^cosh7r_. (2.16) 



3. IR renormalons in Wilson loops 

Expression (2.13) for the scattering amplitude was found by summing all large Sudakov log- 
arithms which were artificially extracted from the uniquely defined scattering amplitude T^y 
(2.1). As a result the perturbative expansion for the scattering amplitude, as we will show in 
this section, is not well defined, and has ambiguities associated with IR renormalons. To restore 
the uniqueness of the physical quantity, the perturbative expression should be supplemented 
by nonperturbative corrections. However, we have to pay for this by introducing a new scale, 
which characterizes the size of nonperturbative effects. We do not know how to evaluate the 
nonperturbative corrections but we may predict in general their structure and dependence on 
impact parameter b by exploiting the idea of IR renormalons. 

The occurrence of infrared renormalons in the amplitude may be seen explicitly in an "im- 
proved" calculation of the Wilson line (WL) expectation value by replacing by the running 
coupling constant as{k'^) [15]. As an example, we consider the one-loop calculation of Wi of 
Fig. 2 in Feynman gauge, and find the structure of the nonperturbative correction for this 
particular diagram. 

W, = 4(i« nil'-''{pm) I 7^ da d/3e^(fi/5+^)^e-(^'^^)"^^^ , (3.1) 

J (^vr) J —oo J —CO 

where D = 4 — 2e , d^k = dk^dk^d^^'^k . After integration over a,f3 and k± in (3.1) we get 
the vacuum average of the WL as an integral over the gluon transverse momenta: 

W, = 4(t« i«)^^-^(-i7rcoth7) / .f'l\rf'^ (3.2) 
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Notice that after integration over k"^, even for fixed coupling constant this expression contains 
infrared poles in e. To regularize IR divergences we have introduced dimensional regularization, 
with scale /i. The IR cutoff A (see Sect. 2) is neglected for this argument, and we may assume 
for a moment that A ^ A. The cross divergence does not appear for nonzero impact parameter 
b, because b regularizes the gluon propagator at short distances. Let us now review how IR 
renormalons appear [15]. 

After substitution of the relation ^^(P) = /o°° da{k'^ / Is?)~'^^'^ into (3.2), and after integration 
over transverse momenta, we get 

M/. = -(f»n-(-ooth,)(./6^)-/j— ^^^^-^(^— j . (3.3) 

We notice that {^^^^'^ ^ e""^/"" with ag = as{'i/b'^) ■ We thus identify the right-hand side of 
(3.3) as the Borel representation 

n(as) = r da7r(a)e-''/''' (3.4) 
Jo 

of 7r{as) = Wi , with 

n{a,) = -{t'^^nk-incoth^){nfl%y [ /wT l/f^ ■ (3-5) 

TT J [e + pia) i(i + piaj 

The limit o" ^ produces a singularity in e due to IR divergences in (3.1), because we have 
neglected the IR cutoff A. Away from o" = we put e = and find that the function tc{(t) 
has singularities generated by F -function at a* = l//3i, 2//3i, 3//3i, . . . . These are the infrared 
renormalons. The result of integration in (3.4) depends on the regularization prescription. 
That is, the factor e"*^*/"" = (6^A^)" induces an ambiguity in the scattering amplitude at the 
level of the power corrections 0(6^A^)". This fact, together with the observation that the 
infrared renormalons come from region of small gluon momenta (A; ~ A) where the coupling 
constant becomes large and where we expect the appearance of nonperturbative corrections, 
implies that the nonperturbative corrections are power corrections. The first renormalon gives a 
contribution at the level 0(b^A^). Thus, for the WL to be well defined, nonperturbative effects 
should contribute at the same level. Hence the nonperturbative correction has the following 
form: 

p(r ®r)(i7rcoth7)(62A2), (3.6) 

where p is some parameter characterizing their size. 

Let us consider the color matrix structure of nonperturbative corrections. The direct product 
of the gauge group generators can be decomposed into the sum of invariant tensors 



Therefore the structure {t"' ® t")(i7r coth7), which corresponds to the particular configuration 
Wi of Fig. 2, reproduces the first two elements of the matrix cross anomalous dimensions (2.7). 
Notice that the structure m coth 7 containing the dependence on 7 is factorized after integrating 
over a and Wilson line parameters and over k±. Evidently the angle and group structure of 
the one-loop diagrams of W2 reproduce the other two elements of the matrix cross anomalous 
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dimension (2.7). Therefore we conclude that the factorized matrix structure of nonperturbative 
corrections is exactly rcross(7) ■ 

In summary, we may represent Wi as a sum of perturbative and nonperturbative terms as 

iy„(7, x%') = (s^ + r:Ls(7) r ^ V - ^^rossh)^M w,{^, i) , (3.7) 

where a = 1, 6 = 1, 2 and A is IR cutoff. We must emphasize that we choose A > A to be able 
to consider the influence of nonperturbative effects. In general the nonperturbative parameter 
p depends on the IR cutoff A. We recall that in this section and afterward A is fundamental 
QCD scale. 

4. Full scattering amplitude 

It is now natural to generalize the exponentiated expressions (2.13), (2.15), (2.16) to include the 
nonperturbative corrections. The Wilson lines (2.2) are defined beyond perturbation theory. In 
perturbation theory, however, we know the renormalization group equation (2.4) for the Wilson 
line and its solution. That is why it is natural to consider the exponentiation of nonpertur- 
bative corrections as well as perturbative. In fact at one-loop level, the matrix structures of 
perturbative and nonperturbative corrections are the same. The origin of this property is the 
following. 

As was shown in the previous section, in Feynman diagrams contributing to the Wilson loop 
nonperturbative corrections come from the integration over small transverse momenta. At the 
same time, the 7 - dependence appears when one integrates over small angles between gluon 
momenta and quark momenta . Since integrations over small angles are completely independent 
of the integration over small transverse momenta, the 7- dependence of both perturbative and 
nonperturbative results coincides. In the paper [28] it was shown that higher order corrections 
preserve the asymptotic behavior of the eigenvalues of the matrix cross anomalous dimension. 
Together, these arguments imply that nonperturbative corrections have r_ ~ I/7 and r_|_ ~ 7 
asymptotics. Therefore the generalization can be performed in the same manner as the one-loop 
example before (3.7) 

T±(P) = / d%e-'' hxp f-r± ['^' - T^U%' ] . (4.1) 

J \ Jx TT T TT J 

For small values of b the scattering amplitude gets its main contribution from the perturbative 
region. At large b the nonperturbative corrections become important. 
For invariant amplitudes we find expressions hke (2.17) and (2.18): 

r(°)=sinh7r_ (4.2) 

r(8)^_^cosh7r_ (4.3) 

We expect that the condition ^ T_ is conserved. Moreover, this condition is provided 
explicitly by the structure ^r_|_A^6'^ in the argument of the exponent, because in high energy 
limit r+ — > 00. 
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Now we come to the main problem of our paper : how to evaluate the integral over impact 
parameter b in the expression for TL (4.1). We will use the method of representation of the 
Fourier transformation via a Mellin transformation, which was so elegantly applied to the 
Drell-Yan process [32]. We should mention also the papers [33], which contain some interesting 
calculations based on Mellin transformation technique. Let us briefly formulate this method. 

If we define 

Tm{s) = J ^{A%yf.{b') , (4.4) 
where T-{b'^) is the Fourier transform of ^-(P): 



f'_(62) = 




e*^r_(^2), 



then the Mellin theorem gives us T^{k^) in the following form: 
where 

r„,) = /^,AVrexp(-r_f ^^M^-r_i^Av). ,4.e) 

Therefore we start our detailed calculation of the scattering amplitude with (4.5) and (4.6). 
We consider two cases: frozen and running coupling constants. 



5. Scattering amplitude for a frozen coupling constant 

First we consider the situation when the coupling constant does not run. For this case, one will 
be able to get an exact answer for the scattering amplitude. Let us rewrite expression (4.6) for 
Tm{s) taking into account that ^^(r) = as- We get 

Tm{s) = J ^(AV)^exp (^logAV - ^r.fe^A^) . (5.I) 
It is convenient to denote 

a 



2n ' 

R^ = (^T_A.\ (5.2) 

TT 

From the beginning let us integrate in 6-space, and after that concentrate our efforts on the 
integration over s. Evidently, Tm{s) is proportional the Gamma-function. Substituting the 
expression obtained for Tm{s) in the expression (4.5) for the scattering amplitude, in which it 
is useful to shift the integration variable s — > s — 1 , one finds 
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Note, if we deal with only nonperturbative corrections {a — aj^ -jli^ = 0) the expression (5.3) 
is simplified. Imagining the contour to envelope the right half plane and summing over all 
residues of r(— s) we find 

— * 

This answer coincides, of course, with the result of direct integration over impact parameter h in 
the expression (4.1) for T_, and sheds some light on the occurance of the Gaussian distribution 
on transferred momentum. 

Let us continue and recall the definition of the confiuent hypergeometrical function (CHF) 
by means of Mellin integral representation (see Appendix) . Taking into account the remarkable 

property of the CHF (A. 2), and identifying a — OL + \^c — \^x — one obtains 

= (^^Je-^F{-a,l,xni + a). (5.5) 
A more detailed form of this expression, using definition (5.2), is 

T iP) = ^(^_!LX^ ,-^W-.f f-^ 1 V(l + ^) (5 6) 

^^"^ pA^r.UvrJ 27r ' ' 4AVr-J ^ 27r ^- ^^-^^ 

This is the exact expression for the scattering amplitude in the case when the coupling constant 
does not run. Let us find the asymptotics of (5.5) using the properties of the CHF (A. 3), (A. 4). 
There are two limiting cases to consider: a; -C 1 and a; ^ 1, and we certainly have a critical 
value of transferred momentum kc^t such that 

4A2 TT " ^ iV3 logs/m2 ' ^ ' ' 

which corresponds to a; = ^ = 1. At this value there is a crossover between rapid Gaussian 
decrease over k in the scattering amplitude and slower Regge dependence in t. 

• 1. /c -C /ccrit, which corresponds to x <^ 1. 
Substituting (A. 3) in (5.5) we get 

^-^^'^ = h (^) '"^^ " "^"^ " ^^^^^ + • • ^^-^^ 

where 'il){x) — dlogV {x) / dx . The result is a Gaussian distribution over the transferred momen- 
tum. 

• 2. A; ^ /Ccrit; which corresponds to x ^ 1. 
Substituting (A. 4) in (5.5) we get 

^ ' '_k'^\w?) Til -a) ^ ' 

The result has the standard reggeized form with Regge trajectory a{t) — —j^Nlog^ and 
-t = P. 
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6. The scattering amplitude as an asymptotic series: Run- 
ning coupling constant. 

Let us substitute the expression for the running couphng constant, Q;s(r) = l/(/3ilog^) in 
(4.6) and integrate over r. Then 

Tm{s) = / ^(AV)^exp(-i^logi^ - ^r_b^A^) . (6.1) 
J ^Trpi log TT 

To evaluate a double integral like (4.5), Collins and Soper proposed in [32] a saddle point 
approximation. For the Drell-Yan process this method is applicable because in the high energy 
limit the nonperturbative corrections play no role, and the saddle point lies in the perturbative 
region. In quark-quark scattering, however, the application of this method is doubtful from 
the beginning because the factor r_, which appears in the argument of the exponent in (6.1), 
approaches zero in the high energy limit. Moreover, the position of the saddle point is very 
sensitive to the nonperturbative parameter p. The saddle point may lie in perturbative region 
or out of it. That is why we suggest a method of evaluating the integral based on the sum of 
the contributions from all poles in s. 

It is useful now to introduce new variables, 

r_ 

n 



27r/3i' 



i?' = -^r.A' , (6.2) 

TT 

Using an a-representation for 7 \ yn and integrating over b in (6.1), one finds 

^«W = rMi (ip) + (6.3) 

The integral Tm{s) is not well defined at ioo because the expression 7^ — is divergent 

at IP' = 1/A^. This is a consequence of the divergence in the running coupling constant at 
P — A^. Taking this remark into account, let us go ahead and try to represent the integral 
as an asymptotic series. We also find the condition under which the asymptotic series can be 
approximated by the first terms and show that it is reasonable. 

In the expression (4.5) for T^{k'^) we shift variable s to s — 1 and substitute the result (6.3) 
for Tm{s). We get, 

T_{k'^)^^^^J'^ dsi^ Pr(l + g + a) i^X d'-^da. (6.4) 

^ ' 2m B?V{n) J-ioo yAR^ j r(s + 1) io ^ / 

We notice that this expression is similar to (5.3) except for the integration over a. Using the 
definition of the CHF by means of Mellin integral representation (A.l), taking into account the 
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property (A. 2) of the CHF and identifying a 



a 



+ 1, c = 1, a; = ^ , one obtains 



r{n) 



e 4fl2 



r(l + a)F(-a, 



a 



n-l 



da . 



(6.5) 



As in Section 5 we concentrate our attention on the two hmiting cases. The critical value of 
the transferred momentum is determined as before by (5.7). 

• 1. <C /ccrit, which corresponds to x <^ 1. 

Substituting (A. 3) in (6.5) we get 



e 1^ r(l + a){l - xa + ...){— j a 



n-l 



da 



TT 



where 



e ( Ji - xJ2 + . . .) , 



r(n) 

pioo 

/ r(i + a) 

^0 



J' 



An" 



a"'da . 



(6.6) 



(6.7) 



(6.8) 



One can easily determine corrections by keeping more terms in the expansion of the CHF. 

Consider the integral Ji (6.7) and represent it as an asymptotic series. In this case, one can 
turn the contour from (0, too) to (0, oo) or change variable a to ia. After expanding the Gamma 
function r(l + a) as a series in powers of a and integrating over a we get 



rW(l) T(n + r) 



Tin) 



r=0 



ri 



, , n^(l) , n(n+l)W/(l)+-i/.^(l)) , 

^ ' . r>2~r . _,9 ~r... 



v 



(logf)' 



(6.9) 

We now have the expression (6.9) for Ji in terms of an asymptotic series. The first terms of 
this series yield accurate value for Ji when 



log —I = llog^ 
^ A2 I I ^ 4A2 



> 1 



(6.10) 



Note that in high energy limit n — > and therefore we didn't include n in the condition 
(6.10). J2 differs from Ji only by an additional power of a in the integral. It means that the 
contributions from J2 are suppressed compared to those from Ji. Evidently the same argument 
can be applied to Ji = r(l + a) a"'~'^~^^da, i = 3, 4,... Therefore we conclude that the 

scattering amplitude can be approximated under the condition (6.10) by the first terms of the 
asymptotic series: 



TT 



9' 



I 



e 4i{2 



1 + 



n 



I 



(V'(l) - x + 0(x2)) + 



+ 



(6.11) 



) ) 
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A more detailed form, using definition (6.2) is 

Therefore we have the Gaussian distribution over transferred momentum. Let us investigate in 
more detail the first exponent in (6.12). Using the expression (2.9) for the eigenvalue r_ one 
can rewrite it as 

P 1 
S \ 4A^ TP 



rn? 



(6.13) 



This means that by increasing energy s the bulk of the diffraction peak, which is concentrated 
for 

^ 2A2 7V2 _ 1 I 

7rp log s/m? ' 

becomes narrower. This phenomena is called shrinkage, which we have thus derived from IR 
renormalon analysis. The slope of the soft gluon trajectory is 

1 1 A^^ 

and we can estimate the value of nonperturbative parameter p using the experimental result 
for the slope: a' = 0.25Gev-2 [35]. It turn out that p ~ 27. 

Let us now substitute (6.12) into expression (4.2) for invariant vacuum amplitude and 
evaluate the differential cross section as da/dt = l/s'^\T^'^^\'^. Taking into account that at high 
energy 2sinh7 — >• s/m'^ and factor l/r_ ~ log(s/m^), we find that the differential cross section 
is given by: 

Notice that the exponential factor in the expression (6.15) has a very interesting feature : 
the nonperturbative parameter p has penetrated into the perturbativc expansion. This factor 
should determine the intercept of the soft pomeron. However, based on our analysis we can 
only estimate it qualitatively. All that we can say is that the intercept depends on the nonper- 
turbative parameter p, depends slightly on the energy through r_, and at high energy is close 
to 1. Moreover, the cross section (6.15) cannot grow faster than log^ s at the limit s — > oo. We 
conclude thus that our soft pomeron satisfies the Proissart bound. 

• 2. k ^ kcrit, which corresponds to x ^ 1. 

Substituting (A. 4) in (6.5) we get the following expression for T_{k'^) , 

T_(^) = r r^(i + f f 1 + <i±i^ + . . .) . (6.16) 



A;2 r(n) Jo 7T J \ X 

Expanding r^(l + a) and sinTro; in powers of a and integrating over a we get an asymptotic 
series 

2r_ 1 
T-{k'') = --5-^^- ^exp 
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which is vahd for 

|log^|»l. (6.18) 

The lowest order is independent of p. This means that the nonperturbative corrections can be 
neglected in this region. The formula (6.17) is a generalization of the standard reggeized form. 

The critical value of momentum transfer kcrit (5.7) depends on the energy s. Taking into 
account the formula for the slope (6.14) we find k^^-^^^ = 1 / {a' log s / m^) . Therefore by increasing 
the energy, the value of /Ccrit decreases. We conjecture that this feature together with the change 
in the behavior of the scattering amplitude at /ccrit describes the crossover region (probably the 
dip region) in the pp differential cross section. Since the scale of nonperturbative effects is 
typically 1 Gev (see, for example [34]), we identify fccrit with 1 Gev. Then the condition (6.10) 
is valid. 

To show the consistency between the expressions for the scattering amphtude in case of 
frozen and running coupling constant one can reduce (6.17) to (5.9) and (6.12) to (5.8). by 
taking the limit (3i ^ 0. 

7. Summary and conclusions 

We began by expressing the scattering amplitude in terms of path ordered exponentials, Wilson 
lines, which admit a nonperturbative definition. 

We then showed that the resummation of soft gluon perturbative corrections, leads to am- 
biguities in the perturbative series associated with IR renormalons. 

By examining the IR renormalon structure we predicted the form of nonperturbative cor- 
rections, using the idea that the ambiguity of the perturbative series is compensated by an 
ambiguity in the determination of nonperturbative correction. We thus observed that nonper- 
turbative corrections are power corrections in the impact variable b. We restored the uniqueness 
of the scattering amplitude adding by "hand" the dominant nonperturbative corrections, which 
are parameterized by a new scale. 

We went on to investigate the interplay between both kinds of corrections, and found the 
critical value of the transferred momentum (5.7) at which the behavior of the scattering am- 
plitude is drastically changed. In each region we represented the scattering amplitude as an 
asymptotic series. For k ^ kcru nonperturbative effects play the main role and this region of 
momentum transfer corresponds to the soft pomeron (See (6.12) and (6.15)), which preserves 
the Froissart bound. For k 3> /Ccrit the large perturbative Sudakov corrections are dominant 
(See (6.17)). 

It is naturally to expect that the nonperturbative parameter p can be expressed as an 
expectation value of a nonlocal operator. Indeed, the expression ^W^// (7, &^A^) must be 
related in some manner with parameter p. In [15] this connection was found explicitly for 
the Drell Yan process. However, in case of q-q scattering we have additional color structures 
and mixing of the Wilson lines, and, we leave the discovery of a similar relation to future 
investigation. 

We may generalize our considerations to study high energy elastic quark-antiquark scatter- 
ing. An antiquark can be treated as quark moving backwards in time with opposite velocity. 
All the above results for quark-quark scattering amplitude can be applied, provided that we 
replace the angle 7 between quark velocities by 7 — ivr — 7 in rcross(7)- Then in high energy 
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limit (7^1) the asymptotics of the eingenvalues of the matrix Fcross do not change. Therefore 
the soft pomeron contributes equally to the pp and pp total cross sections. 

In our approach we took into account the contributions from all diagramms describing the 
interaction of quarks with very soft gluons. The basic of BKFL pomeron is harder gluons 
{t ~ m^). We hope that our paper will help to understand how the two pomerons, soft and 
hard, are related to each other. 

Finally, to identify the physical consequences of the resummed nonperturbative corrections, 
we must embed q-q scattering in p-p scattering and conduct an additional investigation. 
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Appendix A. 

Here we review the main properties of the confluent hypergeometrical function (GHF) [36] 
that are used above. The definition of the CHF by means of a Melhn integral representation is 

T,^ ^ 1 r(c) F(-s)F(a + s) _ 

F(a,c,-x) = TTT-T N — -x^ds. (A.l) 

^ ' ' ^ 27ri F(a) Js-ioo F(s + c) ^ ^ 

We use the following properties: 

1. The continuation formula of Ernst Kummer 



F{a,c,-x) ^ e ""Flc- a,c,x) , (A.2) 



2. The asymptotics of CHF 
For x <C 1 



For X ^ 1 
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Figures: 




Fig.l: Integration paths (a) and (b) entering into the definition of the Wilson fines Wi and 

W2, respectively. 




Fig. 2: One- loop diagram contributing to the fine function Wi. Sofid fine represents the 

integration path, dotted lines denote gluons. 
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